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We show that the ARPES spectra for emission from the bonding as weU as the antibonding Fermi 
surface sheet in Bi2Sr2CaCu20g (Bi2212) possess remarkable site selectivity properties in that the 
emission for photon energies less than 25 eV is dominated by p ^ d excitations from just the 0-sites 
in the Cu02 planes. There is little contribution from Cu electrons to the ARPES intensity, even 
though the initial states at the Fermi energy contain an admixture of Cu-d and O-p electrons. We 
analyze the origin of this effect by considering the nature of the associated dipole matrix element in 
detail and find that various possible transition channels (other than p — * d on O-sites) are effectively 
blocked by either the fact that the related radial cross section is small and/or a lack of available final 
states. Our prediction that ARPES can preferentially sample Cu or O states by tuning the photon 
energy suggests novel possibilities for exploiting energy dependent ARPES spectra for probing initial 
state characters in the cuprates. 

PACS numbers: 79.60.bm,71.18.+y,74.72.Hs 
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I. INTRODUCTION 

The photointensity observed in an angle-resolved pho- 
toemission (ARPES) experiment is the result of mod- 
ulation of the spectral density of the initial state via 
the ARPES matrix element. This modulation can be 
quite substantial in the cuprates and depends in general 
strongly on k|| , polarization and energy of the incident 
photons as well as the energy and character of the ini- 
tial stateiSiiiiSiSiSiS Chuang et al.^, for example, have 
recently exploited the large theoretically predicted dif- 
ferences in the ARPES cross sections of bonding and 
antibonding pieces of the Fermi surface (FS) by tuning 
the photon energy to adduce the doping dependence of 
the bilayer splitting in Bi2Sr2CaCu208 (Bi2212). Asen- 
sio et ali4 find that apparently different looking Fermi 
surface maps obtained in Bi2212 over wide area in the 
(ka;, ky) plane at different photon energies largely reflect 
the effect of the ARPES matrix element. In other ap- 
plications, such as a recent proposal"'^^ to determine the 
spectral function of the boson mediating the Cooper pairs 
in the cuprates, it will be necessary to develop strategies 
for minimizing the effect of the ARPES matrix element in 
order to obtain robust physical results through ARPES 
spectra. It is clear that a good understanding of the na- 
ture of the ARPES matrix element is of considerable im- 
portance in continued development of ARPES as a probe 
of electronic structure of complex materials. 

With this motivation, we examine in this article the 
nature of emissions from the bonding as well as the anti- 
bonding parts of the FS of Bi2212. In considering contri- 
butions to photointensity arising from different angular 
momentum channels and various atomic sites in the unit 
cell, we find that the ARPES matrix element possesses 
remarkable site selectivity properties in that the ARPES 
intensity throughout the (ka;,ky) plane in the 5-25 eV 
energy range is dominated by only the O p ^ d tran- 



sitions. The contribution from Cu sites is quite small, 
even though the electronic states at the Fermi energy in 
Bi2212 contain substantial Cu character. In order to as- 
certain the origin of this effect, the dipole transition ma- 
trix element between the relevant initial and final states is 
analyzed in detail and the factors responsible for this be- 
havior are identifiediSi Our prediction that ARPES can 
preferentially sample Cu or O states suggests novel pos- 
sibilities for exploiting energy dependent ARPES spectra 
to gain insight into character of initial states in complex 
materials. 

The paper is ordered as follows. In Section II we 
briefly discuss the underlying formalism. Section III 
presents and analyzes the nature of the ARPES spec- 
tra. In Subsection IIIA, a typical ARPES spectrum for 
emission from the Fermi energy in Bi2212 is considered 
to orient the reader. Subsection IIIB continues with a 
discussion of site-resolved photointensities over the 5-25 
eV photon energy range for k|| values covering the en- 
tire (k^jliy) plane. The origin of the site-selectivity ef- 
fect is then delineated in Subsection IIIC. Finally, Sec- 
tion IV makes some concluding remarks, including pos- 
sible application of site-selectivity property in gaining 
insight into interesting issues such as the Zhang-Rice 
mechanism^^^"-'^^- in the cuprates and the doping depen- 
dence of the Hubbard U parameter via energy dependent 
ARPES measurementsi^. 



II. OVERVIEW OF FORMALISM 

The methodology used in this study has been discussed 
previously in Ref. |2, to which we refer the reader for de- 
tails; See also Refs. lllll4lll5lll6llr^lTMl9l . Some comment 
in this connection is nevertheless necessary in order to 
meaningfully describe the new results presented here. In 
particular, Ref. shows that substantial insight into the 
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nature of the ARPES photointensitj*2^ resulting from ex- 
citation between specific bulk initial and final states in 
the solid, ijji and tpf, can be obtained in terms of the be- 
havior of the corresponding momentum matrix element, 
< tpf\p\ipi >. We then expand the ip within the KKR 
band structure scheme^fliSi as 

Vi(r) = ^z'cfi?f(r)li(f7), (1) 

where L = {l,m) is a composite angular momentum in- 
dex, ¥^{0.) are real spherical harmonics and (3 denotes 
different basis sites. are expansion coefficients and 

Rf{r) is the radial part of the Bloch wave function on 
site (3. The use of form ^ for the initial and final states 
yieldo^ii8,i9i22 

< i^fim >= E E eJ-'-'C^^:clBi'^,g2,L' ■ (2) 

Q,/3 L,L' 

Here denotes a unit vector along a-direction and 
primed indices refer to the final state. Qf^ j^, are the stan- 
dard Gaunt coefficients. The detailed expression for B^^, 
is given in Ref. i^- For present purposes, the important 
point to recognize is that 

5f;, oc / Rf,{r)* r i?f (r) dr , (3) 

where the integral extends to the muffin-tin radius r^.jx 
of the atom /? in the basis. It is useful to decompose 
the momentum matrix element into contributions from 
various angular momentum channels and sites in the unit 
cell as follows: 

M„ = < I p„ I > = ^ M2 (4) 

7 

= Y.J2Kl,l', (5) 

7 L,L' 

with 

M2, L,L' = E (^L')* ^l' ^M' ^L-L' ■ (6) 

S 

Aside from the obvious notation in Eqs. I4I6I indices 7 
and 5 are together meant to encompass the summation 
of Eq. 121 over all basis sites (3 in the unit cell. In the 
present case of tetragonal Bi2212, 7 takes on eight dis- 
tinct values which include the Ca atom and seven pairs 
of symmetrically located atoms (Bi, Obi, Sr, Osr, Cu, 
Ocu.x, Ocu.y)- S in Eq.|Bl on the other hand, takes only 
two values, which account for the two atoms placed sym- 
metrically with respect to the Ca layer in the unit cell. 
In this sense, we may think of 7 as a "site" index and S as 
a "pairing" index. In this work we study particularly the 
matrix element M2 ^ ^/ , which includes the contribution 
of pairs of sites S related by mirror symmetry. We will 
see that its dependence on different sites and excitation 
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FIG. 1: Grey scale plot (white denotes high) of theoretical 
ARPES intensity for emission from the Fermi energy at hv = 
21 eV for tetragonal Bi2212. Light is assumed to be polarized 
along [110] (white arrow). White and black stars give the ky- 
points (corresponding to different angles P) where the energy 
dependence of the dipole matrix element for the bonding (B) 
and antibonding (A) Fermi surface sheets shown in Figs. |21 
and 121 has been computed. 

channels L — > L' is important for delineating the nature 
of photointensity in Bi2212. 

Concerning relevant computational details, the crys- 
tal potential used is the same as that employed in our 
previous studies of Bi2212 on a body centered tetrago- 
nal lattice and involves 30 atoms per unit cell; See, e.g., 
Refs. IHandjl For definiteness, incident light is assumed 
polarized along the [110]-direction. The damping of final 
states is included via an imaginary part of the final state 
self-energy, = 2 eV. In order to understand the ky and 
energy dependencies of the dipole matrix element, calcu- 
lations were carried out over the photon energy range 
of 5-25 eV for each of the 20 k|| -points considered in 
the (k^;, ky)-plane. This set of k||-points covers both the 
bonding and the antibonding parts of the bilayer split 
Fermi surface of Bi2212. 



III. RESULTS 

A. ARPES spectrum of Bi2212 

For orientation. Fig. ^ shows a typical computed one- 
step ARPES spectrum for emission from the Ep at 21 eV. 
The Umklapp and shadow features^ii are not included in 
this simulation for simplicity. The imprints of the bilayer- 
split bonding (B) and antibonding (A) parts of the FS 
are seen clearly and are in remarkable accord with recent 
ARPES experiments^ii^iLS^. The bonding band gives 
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FIG. 2: Intensity I (Eq.|7|l for k||-points (given by the angle 
f} of Fig. lying along the bonding and antibonding Fermi 
surface sheets in tetragonal Bi2212 over the photon energy 
range of /li^ = 5 — 25 eV. Intensities have been smoothed to 
reflect final state damping effect. Results in (a) and (b) are 
not normalized in relation to each other. 



rise to the X (Y) centered hole sheet B. The antibond- 
ing band is responsible for the electron-like F-centered 
sheet A with a very narrow neck around the M-point^. 
The associated intensity appears somewhat diffuse due 
in part to the presence of a van Hove singularity in the 
electronic spectrum close to M. Note that the polariza- 
tion vector of the incident light (white arrow) breaks the 
symmetry between the irreducible regions T ~ Y — M 
and T — M — X of the Brillouin zone. For this reason, we 
have characterized the k-points on the A and B FS's in 
these two regions by the indicated angle /?, where /3 > 
refers to the upper half and /3 < to the lower half on 
the right hand side of the figure. The intensity from ei- 
ther the A or the B sheet is zero along the nodal line 
r — Y (/3 = 45°), but this is not the case for the sym- 
metrically placed nodal line T — X for f3 = —45°. This is 
due to the fact that even though the momentum matrix 
elements and My are non-zero in general they add 
(for the present polarization) destructively along T — Y 
but constructively along the F — X line. 

As already noted, insight into the ARPES intensity 
can be obtained by considering the momentum matrix 
element < ipf \ p li/ji > which connects the relevant initial 
and final states and the associated intensity 

/ = \A- < ^f\p\4'^ > f . (7) 

Fig. 13 shows intensity computed by using Eq. [7|for a se- 
ries of k|| -points lying along the bonding and antibond- 
ing Fermi surface sheets characterized by the angle /3 of 
Fig. n It is clear that the dipole matrix element varies 
strongly with photon energy and k|| . The main features 
in the intensity from both the bonding and antibonding 
bands are the two peaks at hv around 11 eV and 18 eV. 
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FIG. 3: Contributions to the photointensities of Fig. |5| from 
Cu and O sites in the Cu02 planes in the crystal lattice, (a) 
and (c) give contribution from Cu sites for transitions from the 
bonding and the antibonding initial states respectively, (b) 
and (d) similarly refer to the O-sites. See caption to Fig. |21 
for other details. 



B. Contributions to photointensity from different 
sites 

We next consider the extent to which Cu and O related 
initial state electrons are excited in the photoemission 
process. Fig. |31 gives site (7) resolved contributions to 
the photointensities obtained by taking into account only 
the terms related to Cu or Ocu-sites of the Cu02 planes, 
i.e. for 7 = Cu or 7 = Ocu.x+Ocu,y in Eq. ^ where the 
corresponding contributions from the two Cu02 planes 
are added together. Fig. ^a) shows the Cu contribution 
from the bonding initial state as a function of ky and pho- 
ton energy. The intensity is seen to be practically zero 
indicating that, for emission emerging from any part of 
the bonding FS, Cu sites essentially do not contribute to 
the photointensity. In sharp contrast, contribution from 
Oc„-sites in Fig.|2Ib) is substantial. In fact, the remain- 
ing sites in the crystal lattice (not shown) are found to 
give little contribution, which explains why the intensity 
pattern from Ocu-sites in Fig. Olb) is quite similar to 
the corresponding total intensity in Fig. |2Ia). Much of 
the preceding commentary concerning the bonding FS is 
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also seen to be valid for the antibonding FS by looking 
at Figs. |ni[c), |2Id) and[2Ib) and need not be repeated. 
The results of Fig. 13 thus substantially extend the obser- 
vations made in Ref. 2 at a single k||-value (i.e. M{tt,0) 
point) to show that the ARPES matrix element is domi- 
nated by contribution from Ocu-sitcs in the Cu02 planes 
over the entire FS at photon energies below 25 eV. We 
have carried out limited computations at higher energies 
and find that around = 40 eV both Cu and Ocu con- 
tribute roughly equally to the photointensity. 



C. Why copper electrons do not contribute? 

Eqs. |2] and El make it clear that the value of the mo- 
mentum matrix element depends on the magnitudes and 
phases of several component quantities, namely, the ex- 
pansion coefficients C]^'^ of the initial and final states, the 
Gaunt coefficients Q1 ^/ , and the radial integral terms 

Bj'if . Concerning the character of the initial states of 
interest here, note that the antibonding as well as the 
bonding initial states in the vicinity of the FS in Bi2212 
primarily consist of electrons in the Cu02 bilayers, where 
the Ocu,x, Ocu.y and Cu atomic sites contain approxi- 
mately 80-90% of the weight in the associated Bloch wave 
function^. The primary Cu orbital involved is dr^2_y2^ 
while the most relevant Ocu orbitals are Px and Py. The 
distribution of electrons between different sites and or- 
bitals changes as a function of k|| , but in general the 
electrons &t E k, Ep are strongly concentrated in the 
vicinity of the Cu02 planes. 

Since the initial states of interest are mainly Cu d and 
Ocu P type, from various possible L ^ L' transition 
channels in Eq. (HJ the dipole selection rule A/ = ±1 lim- 
its the L' for the final states to the Cu p and / states 
and Ocu s and d states, respectively. A handle on the 
availability of appropriate final states (i.e. a substantial 
final state coefficient C]^', ) is provided by the partial den- 
sities of states (PDOS's) of the final states resolved into 
different T-channels and sites, shown in panels (a)-(d) of 
Fig. ^ In the case of Cu, Figs.^a) and^b) show that 
PDOS is substantial only in the p channel, the value in 
the / channel being quite small, so that significant con- 
tribution to the matrix element is to be expected only 
from the d — > p transitions. Similarly, Figs.Qfc) and^fd) 
show that for the Ocu-rclatcd final states, the d PDOS is 
substantial so that the p d transitions from Ocu-sites 
will dominate. In short, the lack of Cu /-type and Ocu 
s-type final states in effect "blocks" the Cu d —> f and 
Ocu p — > s channels in the allowed dipole transitions 
such that only Cu d ^ p and Ocu p ^ d excitations 
remain to be of interest in further analysis. 

We next consider the term B'j'lf of Eq. which is de- 
fined by the radial integral of Eq. This term may be 
thought of as the radial cross section associated with site 
(7,(5) for the transition from I to I' orbital. In Fig.^e) 
these radial integrals are presented for the relevant Cu 
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FIG. 4: Left-hand side: PDOS's for final states are presented 
as a function of energy above Fermi energy (Ef) for the most 
relevant angular momentum channels, (a) and (b) show the 
PDOSs of Cu p and /, and (c) and (d) show the O s and d 
final states. Right-hand side: (e) gives the radial integral term 
Bj'^j, of Eq.|H]for indicated I I' transitions from the bonding 
initial state at f3 — 0° (near i\f-point) as a function of final 
state energy. The data with crosses (squares) is for Cu (Ocu) 
sites of the Cu02 planes. In (f) and (g), typical initial (solid) 
and final (dash-dotted) state radial wave functions Rf, (r) and 
Rf{r), respectively, are shown as a function of radius r for 
Ocu p ^ d and Cu d p transitions, respectively. The 
dotted line shows the behavior of the radial overlap integral 
(Eq.OJ. 



and Ocu related transitions. All curves are seen to be 
smooth, slowly varying functions of energy, although ab- 
solute values vary greatly between different transitions. 
The Cu d f and Ocu p d excitation channels are the 
strongest, whereas the Cu d ^ p and OcuP ^ s integrals 
are smaller and will attenuate the related excitations. In 
view of the discussion of the preceding paragraph, weak- 
ness of the Cu d — > p radial cross section is particularly 
important, since the Cu d — >• p channel is the only Cu 
transition channel that could have contributed substan- 
tially to the photointensity. It follows then that the Ocu 
p d is the main transition channel in the formation of 
the momentum matrix element for emission from the FS 
in Bi2212 at low photon energiesj^ 

Insight into the relative magnitudes of the Ocu p ^ d 
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and Cu d p related radial cross sections can be gained 
by considering the behavior of initial and final state ra- 
dial wave functions as shown in panels (f) and (g) of 
Fig. 01 In the case of Ocu p ^ d transition, initial and 
final states in Fig.^Jf) are both nodeless and thus the in- 
tegral Bl'l, does not change sign. By contrast, the wave 
function for the Cu p final state in Fig.^Jg) (dot-dashed) 
possesses two nodes so that when it is combined with 
the nodeless Cu d initial state radial wave function the 
integral oscillates in sign yielding a small value of the 
radial cross section for this transition. This mechanism 
where the nodes in the final state wave function play an 
important role is reminiscent of the arguments made in 
connection with the origin of the Cooper minima2li2Si2& 
observed in photoionization cross sections of atomic sys- 
tems at higher photon energies Putting it all together, 
the discussion of this section shows that even though the 
initial states at the FS in Bi2212 possess considerable Cu 
d character, these states contribute little to photointen- 
sity for hiy < 25 eV due to a combination of two effects: 
The d to f channel is suppressed due to the lack of avail- 
able /-type final states, while the d to p channel possesses 
a weak radial cross sectioui^ 



IV. SUMMARY AND CONCLUSIONS 

We consider the nature of ARPES spectra from Bi2212 
for emission from the Fermi energy. Over the 5-25 eV 
photon energy range, the spectra are shown to display 
remarkable spatial selectivity properties in the sense that 
emissions are dominated by excitation from just the O- 
sites in the Cu02 planes, even though the initial state 
wave functions involved possess substantial Cu charac- 
ter. This selectivity applies to the bonding as well as the 
antibonding FS sheet and holds throughout the (ka;,kj,) 
plane. Insight into the origin of this effect is obtained 
by considering the dipole matrix element for transitions 
between the relevant initial and final states and by de- 
composing this matrix element into contributions arising 
from individual sites in the unit cell in various angular 
momentum channels. Our analysis reveals that, of the 
two possible channels for exciting Cu c?-electrons {d ^ f 
and d p), the d ^ f channel is suppressed by lack 
of available /-type final states, while the d ^ p channel 
is effectively blocked due to a small radial cross section 



(i.e. the term -B^i^ Eq. O. Similarly, of the two pos- 
sible channels for O p-electrons (p — > s and p ^ d), the 
p s channel is suppressed since few s-type final states 
are available, leaving significant contribution to the spec- 
trum from only the O p — > d transitions. It turns out that 
the presence of nodes in the final state wave function 
tends to reduce the radial cross section, which is reminis- 
cent of Cooper minimum type effect in photoionization 
cross section of atomic systems. Limited computations at 
higher photon energies show that the contribution from 
Cu sites increases at higher energies, and that by around 
40 eV, the Cu and O states are excited roughly equally 
in the spectrum. 

Our prediction that ARPES can preferentially sam- 
ple Cu or O states by tuning the photon energy sug- 
gests novel possibilities for exploiting energy dependent 
ARPES spectra for probing initial state characters. An 
exciting example in the cuprates would be to test the 
Zhang- Rice mechanismiiiiS and potentially to deduce ex- 
perimentally the value of the Hubbard U parameter as a 
function of doping. Zhang and Rice noted that large U 
greatly restricts double occupancy of Cu orbitals, so that 
the first holes doped at near half filling - equivalently, 
the states at the top (high energy part) of the lower Hub- 
bard band - should have strong O character. Thus, even 
though ARPES cannot see the upper Hubbard band (and 
hence cannot measure U directly), it should be possible 
to adduce U by measuring the relative Cu/0 character 
of states along the lower Hubbard band, and comparing 
this to predictions of appropriate model computations. 
A direct determination of the doping dependence of U 
would be an important confirmation of recent results and 
predictionsi^. 
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